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Abstract 

A trapping mechanism is observed and proposed as the origin of the anomalous behavior recently 
discovered in transport properties of overdamped ratchets subject to external oscillatory drive in 
the presence of quenched noise. In particular, this mechanism is shown to appear whenever the 
quenched disorder strength is greater than a threshold value. The minimum disorder strength 
required for the existence of traps is determined by studying the trap structure in a disorder 
configuration space. An approximation to the trapping probability density function in a disordered 
region of finite length included in an otherwise perfect ratchet lattice is obtained. The mean velocity 
of the particles and the diffusion coefficient are found to have a non-monotonic dependence on the 
quenched noise strength due to the presence of the traps. 

PACS numbers: 05.60.Cd, 05.45.Ac, 87.15.Aa, 87.15.Vv 
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I. INTRODUCTION 



The existence of chaotic behavior, which is the seemingly random complex motion ob- 
served in deterministic nonlinear systems is now well established. In particular, many ap- 
proaches have been developed for characterizing and understanding the nature of chaotic 
motion 

In addition to chaotic behavior, it has also been shown that deterministic systems can 
exhibit anomalous transport and strange kinetics {2], 3), Q, [5]. In analogy with stochastic 
processes, in the case of normal diffusion , the mean square displacement < x 2 > is 
proportional to time t (< x 2 >~ t), while in the case of strange kinetics 3, 4, f|, < x 2 >~ t 7 , 
with 7 > 2 for enhanced diffusion and 1 < 7 < 2 for dispersive motion. The mean square 
displacement can also have a logarithmic dependence on time, corresponding to 7 = 
10] • Strange kinetics as well as diffusive motion have been observed in both deterministic 



nonlinear systems HQ 



3, Q. 



15] as well as thermal ratchets 161 ] . 



In this paper we concentrate on the dynamics of a deterministic thermal ratchet in the 



presence of a driving force. It has recently been shown [16] that quenched disorder induces a 
normal diffusive kinetics in addition to the drift due to the external drive. Moreover this dif- 
fusive motion is enhanced by higher values of the quenched disorder. If the quenched disorder 
has long-range spatial correlations, diffusion becomes anomalous, and both the correlation 
degree and the amount of quenched disorder can enhance the anomalous diffusive transport 
171 ] . Anomalous transport has been found recently in overdamped systems 18|, ll9|]. In Lin- 
der et al [3] an anomalous coherence is reported and P. Reimann et al [jij find divergence 
on the diffusion coefficient. Although our system differs from those previously reported due 
to the presence of quenched disorder and driving force, the transport anomaly presents some 
similarities that will be discussed below. While anomalous transport in quenched disorder 
ratchets was observed for a range of values of the parameters, the mechanism leading to 
this unusual behavior has not been investigated. Transport properties of ratchets (for a 



recent review of ratchets, see 



Q) i 



is a topic of great current interest, due to the 



possib 



application of these mode 
nanoscale friction 



23 



24 



s for understanding such systems as molecular motors 



21 



251 ]. surface smoothening [26], coupled 



as well as mass separation and trapping at the microscale 28 



Josep 



29. 



ison junctions [271 ] . 



30 



22 



31] . The fluctua- 



tions that produce the net transport are usually associated with noise, but they may arise 
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also in absence of noise, with additive forcing, in overdamped deterministic 



33, 



34 



35 



systems 



32|, 



36]. 



overdamped quenched systems [16J] and in underdamped ratchets 

The aim of this paper is to show that a trapping mechanism is responsible for the observed 
dispersive anomalous transport in an overdamped ratchet subject to an external oscillatory 
drive 37j]. In particular, we show that this mechanism appears when the quenched disorder 



strength is greater than a threshold value. The minimum disorder strength required for the 
existence of traps is determined by studying the trap structure in a disorder configuration 
space. An approximation to the trapping probability density function in a disordered region 
of finite length included in an otherwise perfect ratchet lattice is obtained. We show that due 
to this trapping mechanism, the mean velocity of the particles and the diffusion coefficient 
have a non-monotonic dependence on the quenched noise strength. 

The outline of the paper is as follows: in section II we present the single particle model, in 
section III we define the ensemble and the cumulants, in section IV we present the trapping 
mechanism and discuss its consequences in section V. Conclusions are presented in section 
VI. 



II. MODEL 



The motion of a single particle in an overdamped disordered media is modeled by an 
overdamped ratchet subject to an external oscillatory drive in the presence of a quenched 
noise, using the dynamical equation: 



72 = R(x) + F(t) + G(x) 



(1) 



where, 7 is the damping coefficient, R(x) = —dU/dx is the ratchet force, F(t) is the time 
dependent external force and G(x) is the quenched disorder force. 

The periodic, asymmetric, ratchet potential is modeled by the equation: 



U(x) 



sin(x) — — sin(2x). 



(2) 



with the spatial period A = 2tt, as in previous works [161 . l33l . |34| . l35l . l36l . |38j]. The external 
oscillatory force is given by: 

F(t) = Tsm{ut), (3) 
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where T and u are the amplitude and the frequency of the oscillations, respectively. The 
effects of the substrate randomness is modeled by a quenched disorder term of the form: 

oo 

G(x) =aJ2 t(^l H ( x ~ iA ) - H ( x - (i + 1)A)], (4) 

i=— oo 

where the coefficient a > is the quenched disorder strength, H is the Heaviside function 
and are independent, uniformly distributed random numbers in [—1, 1]. The extension 
to correlated disorder is straightforward. The force G(x) is a piecewise constant force for 
every period of the ratchet potential and gives a reasonably realistic representation of the 
effects of the substrate. 

In order to carry out a numerical solution of Eq. [1] we have carried out a fourth order 



fixed step Runge-Kutta method [39|. Since we are interested in the influence of T and a on 



the transport properties of the system, the remainin g p arameters are set to the following 



values, which were used in previous works [16|, |36|, |37|, |40| : 



7 = = 0.5;lu = 0.1. (5) 

We have carried out numerical solutions of the evolution equation using the following di- 
mensionless variables: 

• The dimensionless position x = x/X, which gives the position of the particle along the 
valleys of the ratchet potential. 

• The dimensionless velocity v = v/v w , with v w = X/T. The mean value of v gives the 
transport velocity of a particle along the ratchet. 

• The discrete sequences obtained by sampling x and v with a sampling period T sa = 
T = 2tv/uj: 

x sa = x(kT) ; v sa = v(kT) , (6) 

with (k=0,l,2,...). Using these variables it is possible to detect synchronization with 
the external driving force. 

A typical trajectory of Eq. [1] consists of an oscillation superimposed on a directed trans- 
port motion with average speed < v >. The particular case of < v >= indicates no 
transport along the ratchet. 
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In a perfect lattice (i.e. a = 0) massless particles remain synchronized over the entire 
[r, u] parameter space 4ol |. The bifurcation diagrams of v sa and < v > as a function of T 
are shown in Fig {TJ for T in the range [0,2]. FigJT] shows that v sa is a monotonic increasing 
function of T and < v > is a stepped function with jumps at specific T values. The meaning 
of these jumps may be understood by considering how the particle's position x varies as a 
function of time t. When T is below 0.96 the particle starts in a potential valley and oscillates 
inside the valley in synchrony with the external driving force, returning every T to the same 
position inside the valley. The particle also has the same velocity. This synchronism explains 
why only one value of v sa is obtained: at every sampling time the velocity of the particle 
has the same value within its oscillatory motion. 

Over the region T e [0.96,1.22], every V value has only one value for v sa but now < 
v >— 1, showing that the particle remains synchronized with the external driving force but 
now it advances one spatial period (one valley) during T. As T further increases the particle 
advances 2 valleys, and then 3 valleys during each T (see the labels (+2, —0) and (+3, —0) in 
FigH]), giving < v >= 2 and < v >= 3, respectively. Furthermore it remains synchronized 
giving only one value of v sa . If T further increases, < v > jumps to a lower value, because 
during the positive half cycle the particle goes forward, crossing several valleys, but it returns 
to one or more valleys during the negative half cycle. This explains the labels in Fig. [IJ 
The motion of the particle remains synchronized with the external force through the entire 
range [0, 2] as it is shown by the single value of v sa . 



III. COLLECTIVE MOTION 



We have studied the evolution of an ensemble of noninteracting particles, uniformly dis- 
tributed over one whole potential valley. The initial position of the particles is given by the 
particle density function: 

p(x, 0) = [H(x - X min ) - H(x - X max )], (7) 

with X max %min T 

The ensemble was allowed to evolve up to a time t = 1000, while the positions of the 
particles were obtained at times = 10k (k = 0, 1, ...) and stored for further analysis. In 
order to perform averages over the realizations of disorder a different quenched disorder 
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sequence was used for each trajectory . In this way, the average over the trajectories also 
includes an average over different realizations of the disorder. 

To characterize the evolution of the packet, the first two cumulants, 

C x = (x k ) ; C 2 = (x 2 k ) - C\- (8) 

and their temporal derivatives, 

(V) = limt -> oo (dCi/dtj ; D = limt -> oo (dC 2 /dt) (9) 

were evaluated at the sampling times as a function of time. Here (V) is the mean velocity 
and D is the diffusion coefficient. In all cases considered, it was verified that all higher-order 
cumulants increase slower than t n / 2 , ensuring that p(x,t) is asymptotically a Gaussian and 
can be determined using the first two moments only. 



IV. THE TRAPPING MECHANISM 



The superposition of the force F(t), with zero temporal mean, and the force R(x) with a 
zero spatial mean value, allows the particles to move at different speeds along the potential. 
Consequently, in spite of the zero spatial mean value of the force R(x) produced by the 
ratchet potential: 

(Rh=^ J R(x)dx = ^[U(x + N)-U(x)] = 0, (10) 

the time-averaged mean value felt by the particles is not zero and is given by: 

ti+i 

(R) T =1 J R[x(t)]dt. (11) 

ti 

This is in fact the reason a sinusoidal driving force produces a positive drift motion when it 
is combined with the ratchet potential. 

A significant consequence of the quenched disorder is the appearance of a trapping mech- 
anism, which arises only in the disordered case (a ^ 0). Traps are a small number of 
contiguous valleys with a negative time-averaged quenched-disorder mean-value, that ex- 
actly compensate the positive time-averaged mean- value of the ratchet potential force. This 
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trapping can be predicted from the synchronization analysis of the perfect lattice case. As 
an example let us consider the case T = 1.65. This T corresponds to the synchronization 
zone (+3,-1) in Fig{TJ When a ^ a particle feels a force that is a combination of the 
ratchet force plus a sinusoidal force with variable amplitude T eq between T — a and T + a. 
Thus, disorder enables the particle reach different zones, as can be seen in Fig. [TJ For 
example, for a = 0.1, the available regions for a particle are (+3, —0) and (+3, —1). Then 
the possible values of < v > are a result of the combination of the positive and the negative 
terms: +3 — = 3 and +3 — 1 = 2. Thus, v is bounded between 2 and 3. For a = a z ~ 0.175 
the available regions become (+2,-0), (+3,-0), (+3,-1), (+4,-1), and (+4,-2). Then 
the possible values of < v > are: +4 - = 4, +4 - 1 = 3, +4 - 2 = 2, +3 - = 3, 
+3 — 1 = 2, +3 — 2 = 1, +2 — 1 = 1, +2 — 2 = 0. Since zero is a possible value, then the 
particle can be localized or trapped. This corresponds to a particle going forward 2 valleys 
during the positive half cycle and going backwards 2 valleys during the negative half cycle. 
Consequently, the trapped particle oscillates inside three valleys in synchrony with the ex- 
ternal driving force. While this analysis is not exact, it provides a reasonable explanation 
for both the minimal disorder strength and the corresponding length K in which particles 
can be trapped. 

In order to determine the trapping probability, we will first define the quenched disorder 
forces G(x) of Eq. [TJ in one of K consecutive valleys, as the coordinate of a i^-dimensional 
disorder configuration space. Possible combinations of disorder are studied in this space 
and each combination is classified either as a trap or a non-trap. For example for V = 1.65 
K = 3; then all possible combinations of three consecutive valleys were studied. The results 
are shown in a K = 3 disorder configuration space in Fig. [2j Note that the trapping region 
in this space changes shape as a increases from 0.2 in Fig. [2(a) to 0.3 in Fig. [2j(c). For low 
values of a there are no traps at all. At a critical value of a = (3k — @3 the first trap appears 
(see FigEk). Let us call this trap the basic trap as it consists of K consecutive valleys with 
equal disorder strength (for the case r = 1, 65, K — 3 and (3^ = 0.136). 

Note that the volume of the trapping configuration region further increases with increasing 
a (see Fig. [2b). For a = Pk-i = P2 — 0.21 (see Fig. [2fc), two arms appear, corresponding 
to traps of only K — 1 consecutive valleys (for the case V = 1.65, K — 1 = 2). If a is further 
increased, there exists a higher value (3k-2 = Pi over which traps of only one valley appear. 
Note that the minimum a for which traps appear is (3k 
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When K = 2, the disorder configuration space is two-dimensional. Traps appear when 
a > Pk = @2- in the interval Pi < a < Pi, the probability space (that corresponds to 
trapping events) gradually mutates from a single triangular shape into a square as a grows. 
When a > Pi, two narrow, rectangle-shaped arms appear, because it is possible now for a 
particle to get trapped in just one valley. We note that, in Fig. [2J in which a K = 3 case 
is considered, the Pi value is outside the plotting range. If Pi had been included in the plot 
three rectangular prisms would have appeared, one for each random variable. For instance, 
the prism that corresponds to valley number 1 would have occupied the volume: 



• -a < a£i < -Pi » 

—a < a^2 < a 
^ —a < a^2 < a j 



(12) 



For each a value the cumulative probability of trapping in a .fT-length trap, called p(ot, K), 
is evaluated as the ratio between the volume of the trapping region and (2a) K . 

When a disordered region of finite length L > K is considered the cumulative probability 
of trapping is approximately given by: 

P{a,K,L) = l-Q(a,K,L), 
Q(a,K,L) = {l-p(a,K)) L - K+1 . (13) 

The approximation used to obtain Eq. [12] is based on the assumption that the events of 
actual trapping of the particle in the neighborhood of the K- length traps are independent. 

Let F(a, L) be the fraction of particles traversing the length L of the disordered region. 
We define f(a,L) by the relation: 

f(a,L)=d[l-F(a,L)]/da. (14) 

In Fig. [31 f{<y, L) is compared with dP(a, K, L)/da with the values of P(a, K, L) obtained 
from Eq. [13j The good agreement between the two results confirms that the independent 
events approximation used in Eq. [T3]is clearly valid for a values below Pk-x- Note that the 
local maximum of both curves in Fig. [3] occurs at Pk-i — 0.21, corresponding to the value 
at which the "arms" begin to appear in the disorder-configuration space (see Fig. [2b). We 
note that the value Pk-2 — 0.5720, where K — 2 valley length traps appear, falls outside the 
range of the plotted values. 
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As the disorder region length L grows, the shape of Figure [3] becomes thinner but its left 
end is still at a = (3k (in this case, a = (3$). If L tends to infinity the probability that a 
particle finds a K- length trap tends to unity when a = (3k- Then, the shape of figure [3] 
becomes a Dirac delta function, which is also predicted by Eqs. [T3] and dH 

We have carried out similar studies for many other values of V in the range [0, 2]. We have 
found that in most of the cases studied, the correlation effects are not important and Eq. [13] 
is accurate enough for determining the probability density function P(a, K, L). There are 
a few T values, however, for which the correlation effects cannot be easily neglected, since 
[3k — Pk-i and the "arms" in the disorder configuration space (like the one in figure [2]) 
appear early on. 

As seen in Fig. 4, Eq. [13] is not accurate when a » (3k- This is because the cor- 
relations cannot be neglected in this region. We propose another method, which we call 
the Conditional- Probability Method or CPM, that takes the correlations into account. This 
method is always accurate provided that a » (3\ but for some V values it also works well 
in the entire a > f3\ range. 

Calculations based on CPM involve the following steps: 1) Approximate the probability 
space volume occupied by trapping events (both black and gray dots in Fig. 3), by a 
polyhedron in which the faces are perpendicular to the coordinate axes. The a > (3\ case 
must be considered and rectangle prisms are the main volume component. The smaller 
volume components are optionally considered, if the method's range of valid a values is to 
be increased. If K = 2 we have a probability space enclosed by a surface with faces that are 
perpendicular to the coordinate axes. 2) Write the volume equations of the non-trapping 
events, that is the volume outside the polyhedron, as the union of the set of volumes. The 
probability of having no traps is proportional to this volume. For example, let us consider 
the T = 1.35 case for which K = 2. The non-trapping probability Q(a, K, L) = Q(a,2,2) 
agrees approximately with: 




(15) 



If the smaller surface of trapping events approximately given by 



-(3 1 < aii < -fl 2 , -(3i < <2 < -#2 



(16) 



is taken into account we can instead write: 
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Q(a,2,2)(4a 2 ) 



-fa < a£i < a 
-fa < a£ 2 < a 



U 



— fa < a£i < a 
-fa < a& < -fa 



U 



-fa < aft < -/3 2 
-/3 2 < at£ 2 < at 



(17) 



3) Increase L, sequentially, in single steps, and find the new volume as the intersection 
of the previous volumes. That is, Q(a,K,L — l).(2a)^ L ~ 1 - 1 , expanded for all a^L and the 
original if-dimensional volume expanded for all a£i, a£ 2 , a^-K values. Continuing 
with this procedure and using the most accurate volume equation, we can write the L = 3 
volume in the following way: 



Q(«,2,3)(8« 3 ) 
/ 



' -fa < a£i < 



-fa < a& < a 
, —a < a£ 3 < a j 



U 



n 



LV 



—a < a£i < a 
-fa < a£ 2 < a 
-fa < a£ 3 < a 



U 



' -fa < a£i < a » 

-A < a£ 2 < -fa 
, —a < < a j 

l 



U 



' —a < a£i < a » 

-/3 2 < a£ 2 < a 



u 



1 fa < a£l < -fa\ 

-fa < a& < a 
, —a < < a j 

\ 



—a < a£i < a 
fa < a£ 2 < -fa 
\ -fa < a£ 3 < a J 



(18) 



Applying distributive property, the intersection of volumes is computed, yielding: 



' -fa < a^i < a » 



Q(a<, 2, 3)(8a 3 



—fa < a£ 2 < a 
-fa < a£ 3 < a 



1 -fa < a£i < a » 



' -/3 2 < a£i < a ^ 



-A < a£ 2 < -fa 
y -fa <a£ 3 <a J 



[j -fa < < 2 < a (J 
y-p 2 <ai; 3 <aj \-fa < a£ 3 < -fa j 

-fa < a£i < -fa 
[J -fa <a&<a [J 
, -fa <a£ 3 <a j 

To find the L = 4 hyper- volume, the surface in Eq. [T7] should be expanded to 4D, all 
over —a < a£\ < a, —a < a£ 2 < a as follows: 



^ (-fa < «Ci < -A^ 



—A < ct£ 2 < a 
\-fa < «6 < -A I 



• (19) 
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/ —a < a£i < a \ 

—a < a£ 2 < a 
-f3 2 < a£ 3 < a 



U 



U 



/ —a < a£i < a \ 
—a < a^2 < a 
-Pi < a£ 3 < -fa 
\ -f3 2 < a£ 4 < a J _ 



(20) 



/ —a < a£i < a \ 
—a < a^2 < a 

—@2 < 0£%3 < a 

\-/3i < < - fa/ 

and the volume in Eq. [TH] should also be expanded to 4D all over —a < < a. Finally 
the intersection between these 4D hyper-volumes must be found. 

As L grows, the number of component hyper-volumes in the set increases and their 
intersections are hard to calculate by hand. We calculated the resultant hyper- volume up to 
L = 15 by using a binary method in which intersections are 'and' boolean operators. Fig. H] 
shows that there is excellent agreement when a > f3\ for L = 15. 



V. CONNECTION BETWEEN TRAPPING AND TRANSPORT PROPERTIES 

The presence of the traps has a significant macroscopic consequence in multiparticle 
systems on the variation of the current (V) and the diffusion coefficient D. In order to 
explain this consequence, consider the cumulants C\ and C*2, shown in Fig. [5J as a function 
of time for a packet of 2200 massless particles with T = 1.65, u — 0.1, 7 = 1, L = 50 and 
(jl = 0.5. The corresponding derivatives < V > and < D >, are plotted as a function of 
a in Figs. [6^ - [6)3. We note that as a function of a, four regions with different transport 
properties may be recognized in these figures: 

1) Very low disorder strength region, a < 0.08. In this region the disorder has no effect 
(see the case a = 0.06 in Fig. ED, the first cummulant C\ is proportional to t and there is 
no diffusion as evidenced by the constant value of C2 and the zero value of < D > (see Fig. 
[6fc). Particles can only reach the region (+3, —1) of Fig. [1] and the only possible value of 
the mean velocity of the particles is 2. Thus, < V > remains constant equal to 2 and the 
dynamics is essentially the same as in the perfect lattice case. 

2) Intermediate disorder region, 0.08 < a < (3k = 0.136 (see the case a = 0.12 in Fig. 
[5] and Fig. [6b, where this region is enlarged). In this region, the first cummulant C\ is 
proportional to t, but now C2 ~ t, indicating the existence of normal diffusion. The slope 
of the packet mean velocity changes abruptly for each value of a, where a new mode is 
reached, as seen in Fig. [TJ At a > 0.08 the region labeled (+3, —0) is reached and the mean 
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velocity may have one of the two values +2 or +3. Consequently V increases with a, as can 
be seen in Fig. [6K At a = 0.1 region (+4, —1) is also available and the mean velocity can 
have any of the values +2, +3 and +4. Then V continues increasing with a, with a higher 
slope. Fig. Eb shows an enlarged view of a region of Fig. [6b, where the normal diffusion can 
be observed. As can be seen, for a < 0.08 transport is not diffusive. At a = 0.08 normal 
diffusion starts and as a is increased further, diffusion is enhanced and the same critical 
value of a appear as in the case of V (Figj6h.). 

3) The region a ~ @k = 0.136 (see, for example, the case a = 0.145 in Figj5j) . In this 
region the trapping mechanism has already started. The cumulants C\ and C% increase as 
t H with H > 1, indicating a super-diffusive behavior. The time-dependence of the second 
moment is very complicated during the trapping process and it strongly depends on both 
the quenched disorder realization and the strength of the disorder. This transitory time 
is considerably reduced as a increases beyond the threshold. The trapping mechanism 
produces an abrupt descent in V. The other values Pk-i = P2 and (3k -2 = A do not appear 
in Figsj6] because the probability that a particle gets trapped in a K = 3 valley length trap 
approaches 1 for a disorder region with L = 50>>i^ = 3, and most particles get trapped 
in K = 3 length valleys, regardless of whether K — 1 = 2 and K — 2 = 1 length traps exist 
or not. In Fig. [6b the super diffusive region is clearly recognized by the high values in D. 
Below the threshold value (3k diffusion is normal and D is much smaller than in the super 
diffusive region. Note that C2 varies with time in a very complicated way as long as there 
exist both, trapped and untrapped particles. Untrapped particles suffer a normal diffusion 
process but trapped particles cause the packet to get wider as it evolves, increasing D. In 
fact, for t — > 00 (and consequently L — > 00) all particles get finally trapped and D — > . 

4) The region a^> {3k (see the case a = 0.165 in Figj5]). In this region, complete trapping 
occurs even for small disordered zones and both < V > and D decrease to zero. The system 
undergoes a transition at a = (3k between two different transport regimes: for a < (3k there 
is normal diffusive transport while for a » (3k both transport and diffusion disappear. In 
the neighborhood of a = (3k anomalous diffusion is present. The anomalous transport found 
recently in an overdamped tilted potential model with thermal noise but without quenched 
disorder and driving force 19j . a similar transition between normal diffusive transport is 
given by the variation of the tilting force F and anomalous diffusion appears near the 
critical value F r . 
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The anomalous transport effect of trapping is robust under thermal fluctuations in the 
sense that thermal fluctuation amplitudes of the order of the quenched disorder strength are 
required to destroy this effect. 

VI. CONCLUSIONS 

A novel trapping mechanism is discovered which is proposed as the origin of the anoma- 
lous transport in a multi-particle overdamped disordered ratchet. It is found that once a 
particle reaches a trap it remains localized inside a small region, oscillating synchronously 
with the external force. By means of the disorder configuration space, critical values for the 
disorder strengths were determined and the fraction of particles traversing a disordered re- 
gion were obtained by means of two methods. In the first method, valid for low disorder, the 
effects of correlations between the contiguous traps were neglected. In the second approach, 
which is valid for disorder strengths over all the critical values, correlations were considered. 
The probability density function shows excellent agreement with the simulation data. The 
trapping mechanism presented here explains the singular behavior of velocity and diffusion 
with disorder in overdamped ratchets reported in [16j. This analysis may be helpful in the 
study of other systems exhibiting strange kinetics and also for practical applications such 
as the design of particle separation techniques in multi-particle systems. 

VII. ACKNOWLEDGMENTS 

This work was partially supported by CONICET (PIP 5569), Universidad Nacional de 
Mar del Plata and ANPCyT (PICT 11-21409 and PICTO 11-495). 

VIII. BIBLIOGRAPHY 

[1] S. H. Strogatz, Nonlinear dynamics and chaos: With applications to Physics, Biology, Chem- 
istry, and Engineering. (Addison Wesley, Reading, MA, 1994). 
[2] H. Scher, M. F. Shlesinger and J. T. Bendler, Phys. Today 44, No. 1, 26 (1991). 



13 



[3] M. F. Shlesinger, G. M. Zaslavsky and J. Klafter, Nature (London) 363, 31 (1993). 

[4] J. Klafter, M. F. Shlesinger, and G. Zumofen, Phys. Today 49, No. 2, 33 (1996). 

[5] A. Lichtenberg and M. Lieverman, Regular and Stochastic Motion. (Springer, New York, 1983). 

[6] T. Geisel and J. Nierwetberg, Phys. Rev. Lett. 48, 7 (1982). 

[7] J. A. Blackburn and N. Gr0nbech-Jensen, Phys. Rev. E 53, 3068 (1996). 

[8] E. Barkai and J. Klafter, Phys. Rev. Lett. 79, 2245 (1997). 

[9] E. Marinari, G. Parisi, D. Ruelle, and P. Windey, Phys. Rev. Lett. 50, 1223 (1983). 

[10] J. Krug and H. T. Dobbs, Phys. Rev. Lett. 76, 4096 (1996). 

[11] R. Morgado, F. A. Olivera, G. G. Batrouni and A. Hansen, Phys. Rev.Lett. 89, 100601/1-4 
(2002). 

[12] K. Mallick and P. Marcq, Phys. Rev. E 66, 041113 (2003). 

[13] E. K. Lenzi, R. S. Mendes and C. Tsallis, Phys. Rev. E 67, 031104 (2003). 

[14] N. Korabel and R. Klages, Phys. Rev. Lett. 89, 21402 (2002). 

[15] H. Kunz, R. Livi and A. Slid, Phys. Rev. E 67, 011102 (2003). 

[16] M. N. Popescu, C. M. Arizmendi, A. L. Salas-Brito and F. Family, Phys. Rev. Lett. 85, 3321 
(2000). 

[17] Lei Gao, Xiaoqin Luo, Shiqun Zhu and Bambi Hu, Phys. Rev. E 67, 062104 (2003). 

[18] B. Linder, M. Kostur and L. Schimansky-Geier, Fluctuation and Noise Letters 1, R25 (2001). 

[19] P. Reimann, C. Van den Broeck, H. Linke, P. Hanggi, J. M. Rubi, and A. Perez-Madrid, Phys. 

Rev. Lett. 87, 010602 (2001). 

[20] P. Reimann, Physics Reports 361, 57 (2002). 

[21] R. D. Astumian, Science 276, 917 (1997). 

[22] R. D. Astumian and P. Hanggi, Phys. Today 55, 33 (2002). 

[23] L. Daikhin and M. Urbakh, Phys. Rev. E 49, 1424 (1994). 

[24] C. Daly and J. Krim, Phys. Rev. Lett. 76, 803 (1996). 

[25] M. R. Sorensen, K. W. Jacobsen and P. Stoltze, Phys. Rev. B 53, 2101 (1996). 

[26] A. L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth. (Cambridge University 
Press, UK, 1995) 

[27] I. Zapata, R. Bartussek, E. Sols and P. Hanggi, Phys.Rev. Lett. 77, 2292 (1996). 

[28] L. Gorre-Talini, J. P. Sspatz and P. Silberzan, Chaos 8, 650 (1998). 

[29] I. Derenyi and R. D. Astumian , Phys. Rev. E 58, 7781 (1998). 

14 



[30] D. Ertas, Phys. Rev. Lett. 80, 1548 (1998). 

[31] T. A. J. Duke and R. H. Austin, Phys. Rev. Lett. 80, 1552 (1998). 

[32] P. Hanggi and R. Bartussek Lecture Notes in Physics Vol. 476, ed. J. Parisi et al, (Springer, 
Berlin, 1996). 

[33] P. Jung, J. G. Kissner and P. Hanggi, Phys. Rev. Lett. 76, 3436 (1996). 

[34] J. L. Mateos, Phys. Rev. Lett. 84, 258 (2000). 

[35] M. Barbi and M. Salerno, Phys. Rev. E 62, 1988 (2000). 

[36] C. M. Arizmendi, F. Family and A. L. Salas-Brito, Phys. Rev. E 63, 061104 (2001). 
[37] F. Family, H. A. Larrondo, D. G. Zarlenga and C. M. Arizmendi, J. Phys.: Condens. Matter 
17, S3719 (2005). 

[38] H. A. Larrondo, F. Family and C. M. Arizmendi, Physica A 303, 78 (2002). 

[39] W. H. Press, S. A. Teikolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C. 

(Cambridge University Press, Cambridge, 1995). 
[40] D. G. Zarlenga, H. A. Larrondo, C. M. Arizmendi and F. Family, Physica A 352, 282 (2005). 



15 



3.5 
3 

2.5 
2 

>* 
h— > 

O 1.5 

0) 

> 

1 

0.5 


-0.5 



~i 1 r 



— V 

sa 








(+3,-0) 


(+4,-1) 


< V > 
























(+2-0) 








(+4,-2) 



(+0,-0) 




(a) 



J L 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

r 



FIG. 1: Sampled velocity v sa , and mean velocity (v) of a particle in a perfect lattice, as a function 
of T. The particle starts at x = 0. Note the jumps in {v} at T ~ 0.96, 1.22, 1.47, 1.57, 1.75, 1.95, 
but i; sa has no bifurcations in this range of T. The label over each zone indicates the number of 
valleys crossed by the particle, forward (+) and backwards (-) in a period T. 
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FIG. 2: Three dimensional disorder configuration space for T = 1.65, u = 0.1, 7 = 1 and fi = 0.5 
as a function of a. The axes represent the quenched disorder in three consecutive valleys. The 
black dots correspond to triplets acting as traps, (a) The case a = 0.20 <; (3k- Note the small 
trapping region at the corner a£\ = a^2 = = ~Pk, (b) The case a = 0.22 ^ (3k-i- Note the 
appearance of two dimensional arms at a = 0k-i\ (c) The case a = 0.30 > @k-i 
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FIG. 3: The probability density function dP(a, K, L)/da from eq. [13] (bold curve) and the function 
f(a, L) from Eq. [TH (thin curve) for 2200 massless particles traversing a L = 50 disordered region, 
are plotted as a function of a, for T = 1.65,7 = 1, w = 0.1 and fi = 0.5. 
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FIG. 4: The fraction of particles traversing a length L = 15 disordered region is plotted as a 
function of a, for T = 1.35, u> = 0.1, 7 = 1 and = 0.5. Circle markers corresponds to simulations 
obtained with 1600 paticles. Dotted line corresponds to function Q(a, K, L) in Eq. [T3l where no 
correlation is considered between contiguous traps. Solid line corresponds to Q(a, K, L) found by 
using CPM. 
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FIG. 5: The cumulants (a) C\ and (b) C2 are plotted as a function of time t, for 2200 massless 
particles, traversing a L = 50 disordered region, with T = 1.65, a = 0.15, uj = 0.1, 7 = 1 and 
H = 0.5. 



20 



3r 
2.5 

2- 
1.5 



(a) 



p=0.?36 



<V> 



1 

0.5 



2000 



T=1 .65 
0.05 



0.1 

a 



0.15 



0.2 




FIG. 6: (a) The mean velocity (V) and (b),(c) the diffusion coefficient D are plotted as a function 
of a (for a packet of 2200 particles traversing a L = 50 disordered region with T = 1.65, lo = 0.1, 
7 = 1 and = 0.5). (c) is an enlargement of part of Fig. (b), where trapping does not exist and 
the transport is diffusive 
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